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I. INTRODUCTION 



In systems of correlated fermions on a lattice some interesting and also puzzling physics seems to happen when 
interaction-induced localization tendencies, antiferromagnetic, and superconducting fluctuations get mixed. The 
standard example of such a system are the copper-oxide superconductors.! In the underdoped regime, between the 
AF and SC phases, correlations of both AF and SC type are strongly enhanced, and a pseudogap is visible in the 
one particle spectrum and in the spin response functions. The pseudogap regime is limited from above by a crossover 
temperature T co (x), a monotonously decreasing function of doping. At the temperatures T > T co the underdoped 
materials are "strange metals" : many physical properties are unlike those of a standard Fermi liquid!. In the overdoped 
regime T co {x) merges with the critical temperature for superconductivity and the regime T > T c is merely a Fermi 
liquid. Another interesting feature of the phase diagram is the unusual form of the order parameter. After a rather 
long period of controversies the d x 2_ y 2 -symmetry is finally generally accepted.! This is one of the reasons to believe 
that the pairing mechanisms are tightly related to the antiferromagnetic tendencies and not to the standard phonon- 
exchange mechanisms. The d x 2_ y 2-{ovva of the superconducting correlations subsists also in the pseudogap regime as 
is seen in recent angle-resolved photoemission! and tunneling! experiments. The simultaneous existence of strong AF 
correlations, as seen by NMFt! or neutron-scatteringi experiments, and even localization tendencies as flattening of 
the band! make us conclude that the interpretation of this regime only terms of superconducting or antiferromagnetic 
fluctuations only is not sufficient, especially because we expect that they are coupled. 

It is striking that some other apparently completely different systems of correlated fermions have very similar 
properties. A phase diagram with the superconducting phase in the vicinity of the spin density wave (i.e. antiferro- 
magnetic) instability characterizes also the quasi-one-dimensional Bechgaard salts! and the quasi-two-dimensional 
organic superconductors of the ET family,0 where instead of doping the relevant parameter for the phase diagram 
is pressure. However the common feature of all these compounds is that they are systems of correlated fermions 
with reduced dimensionality (D < 3) and with strongly anisotropic and more or less nested Fermi surfaces. The 
main points for the understanding of the three groups of compounds are: (i) the destruction of the nesting by doping 
(cuprates) or by applying pressure (Bechgaard salts and ET) ; (ii) the suppression of the Umklapp processes by doping 
the the half-filled band (in the cuprates and ET-s) or by making the half-filled band effectively quarter-filled through 
the breaking the longitudinal dimerization by pressure (in the Bechgaard salts) . Concerning the Bechgaard salts it is 
interesting to remark that some very recent interpretations of the phase diagram of the (TMTSF)PFg materialS sug- 
gest that the intermediate regime between the high-temperature ID behavior and the low temperature 3D physics is 
a strange 2D-liquid with properties very similar to those of the underdoped cuprates above the crossover temperature 

CO- 

From the theoretical point of view it is certainly interesting to construct a theory able to treat on the same 
footing antiferromagnetic and superconducting tendencies in more than one dimension and to follow how the result 
changes with some external parameter that destroys nesting and the Mott-like localization. The first question one 
can ask is whether a purely repulsive model like for example the Hubbard model (or some generalization of it) 
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contains the coexisting and inter-depending antiferromagnetic and superconducting correlations. In such a model 
the antiferromagnetic fluctuations are associated with the enhancement of the particle-hole (p-h) propagators at low 
energies and the superconducting tendencies appear through particle-particle (p-p) propagators. The Hubbard model 
is appropriate because at half filling already a simple mean-field calculation gives an antiferromagnetic instability 
at a finite temperature. However if one tries to include also the p-p processes already in the weak coupling limit 
the problem becomes nontrivial already in the weak coupling limit: a simple mean-field theory is not able to follow 
both p-h and p-p correlation channels. One can of course try to remedy this problem by including summation of 
selected subseries of higher-order diagrams. One such attempt was to calculate the effective Cooper amplitude as 
a sum of bubble and ladder RPA seriesEBE^. The resulting Cooper amplitude is then used as coupling constant for 
assumed effective BCS theory. This procedure thus explicitely decouples three diferent summations (RPA bubbles, 
RPA ladder, and BCS ladder) without real justification. The FLEX (conserving) calculations^ based on the similar 
simplifications are also prejudiced by the choice of diagrams to be summed. The only way to proceed systematically 
is to construct the renormalization group that takes into account all p-p and p-h loops of a given order (or to use the 
equivalent parquet approach). In (quasi-) one dimension the renormalization group has been successfully used and is 
one of the basic theoretical ingredients in the physics of low-dimensional metalsiEl H 

In two dimensions only a limited number of simplified cases was solved by the renormalization group. The poor man's 
scaling applied only to the interactions between electrons placed at the van Hove points it gives an antiferromagnetic 
instability at half- filling and superconductivity of d x 2_ y 2 symmetry if the deviation of the chemical potential /i from 
its value at half-filling becomes of the order of critical temperature of the antiferromagnetic state.0 The equivalent 
parquet approach has been used for half-filling (but without the limitation to the van Hove points) and also finds 
an antiferromagnetic instablity.SB Parquet calculations for simple flat Fermi surfaces@ give an antiferromagnetic 
instability but can not provide a continuous phase diagram as function of some imperfect nesting parameter or band 
filling: the d a .2_ J ,2-like superconducting pole appears simply by cutting the p-h part of the flow, as in ref. [Tt]. The 
scaling approach to a system with a Fermi surface with both flat and curved partsS has also reported a superconducting 
instability in a purely repulsive model, together with deviations from Fermi liquid behavior. However, the complete 
one-loop renormalization group (or parquet) for the real band of electrons with imperfect and tunable nesting, or 
doping, still remains unresolved. The main difficulty is related to the correct treatment of the coupling between p-p 
and p-h channels. 

Different authors have tried to avoid to take into account the coupling between the different renormalization channels 
making drastic simplifications or limiting themselves to some particular forms of the Fermi surfaces or only to the low 
energy effective action. In our former publicationslllsl we have shown that in the Hubbard model one can treat the 
p-h channel perturbatively if the filling is sufficiently far from one half. Then the renormalization group gives only a 
weak Kohn-Luttinger like pairing. The p-p part of the flow is decoupled from the p-h one in the low energy regime 
for the simple reason that the p-h part is negligible there. Other calculations based on the perturbative treatment of 
the p-h channel were also reported.@@ If the Fermi surface is well (but imperfectly) nested, and this is exactly the 
interesting regime, this strategy does not work any more because both p-p and p-h loops are non-perturbatively large, 



3 



even for weak interactions. In the case of the square Fermi surface (with or without the van Hove singularities) taking 
only the leading logarithmic part the coupling between the channels into accountB is equally insufficient. Another 
way to proceed is to see the 2D Hubbard system as an ensemble of coupled chains :0 this approach gives a phase 
diagram with superconductivity formed by pairs of electrons on different chains, giving rise to a spatially anisotropic 
version of an d-wave order parameter. Among the number of the theoretical approaches to the Hubbard model other 
then via the loop-summations the Monte Carlo calculations take in principle "everything" into account but it is still 
unclear whether they giveS or notHi the superconductivity. 

In the present paper we search to reliably determine the phase diagram of the Hubbard model in the vicinity 
of half-filling where p-h processes are non-perturbatively enhanced and at least nearly as important as p-p ones. 
We also detect the dominant components of the angle-resolved correlation functions for antiferromagnetism and 
superconductivity as function of temperature. This allows us to know the symmetry of the microscopic fields whose 
fluctuations become important. The method that we will use is a generalization of Shankar's renormalization-group 
approach0to an arbitrary form of the Fermi surface. In particular, the Kadanoff- Wilson mode elimination (developed 
by Shankar for 2D fermions) applied to the effective action with only two-particles interaction keeps only the strictly 
logarithmic contributions to the flow. Thus even if the nesting is very good but not perfect the p-h part of the 
flow would be zero because the logarithmic singularity is destroyed by imperfect nesting. To keep the p-h part of 
the flow finite even in the case of imperfect nesting we start by formulating the exact Kadanoff- Wilson-Polchinski 
renormalization group for fermions on a lattice. It was formulated previously0El only for the quantum fields with one 
zero-energy point in the momentum space, like the cf> field theory (critical phenomena). In many-fermions system 
in more than one dimension we have on the contrary the whole Fermi surface that plays the role of the zero energy 
manifold, what makes the calculations more complicated. Starting with the full bandwidth as the initial energy cutoff 
we perform iterative mode elimination reducing the cutoff A around the Fermi surface. Collecting at each step of the 
renormalization all the terms (cumulants) of the effective action we obtain the Polchinski equation for the vertices of 
the effective theory at the given step of the renormalization. It is important that even if the initial interaction was 
only a four-point function (two-particle interaction) vertices of all higher orders are created by the renormalization 
procedure. Once the exact renormalization group is formulated we proceed with its truncation at the one-loop level: 
the one-loop truncation of the flow for the four-point vertex is done by neglecting all renormalization-group-created 
vertices of order larger than six. ShankarS already has remarked that the six-point function created by the mode 
elimination is essential to get non-logarithmic contribution to the four point vertex (the effective interaction). The 
one-loop renormalization of the interaction that we obtain in this way appears to be generally non-local in A, i.e. the 
flow of the vertex at a given step of the renormalization depends on the values of the vertex at former steps. This is 
certainly not a pleasant but (as far as we can see) necessary property of the KWP procedure if we want to keep more 
than just the purely logarithmic contributions. 

We than apply the one-loop KWP renormalization group to the Hubbard model. One further approximation we 
make is to consider the effective interaction as a function only of the projection of the momenta to the square Fermi 
surface (marginal interactions) while the radial dependence and dynamics are neglected because they are irrelevant 
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with respect to the Fermi-liquid scaling.ES We neglect also the renormalization of the self-energy. If we take only 
marginal interactions into account this is justified at the one-loop level because the renormalization of the weight 
and of the lifetime of the electrons receives a nonzero contribution only at the two-loop level. We thus renormalize 
only the interaction U{8\ 1 82, O3), a function of three angular variables corresponding to the angular parts of the 
three external momenta, the fourth being determined by momentum conservation. We allow the ^-variables to be 
anywhere on the almost square Fermi surface and not only in the configurations that give perfect nesting or zero 
center-of-the-mass momentum: these two classes of the configurations would correspond only to the processes with 
the leading logarithmic renormalization in the p-h and p-p channels, respectively. As will become clear later, taking 
all three ^-variables without constraints is the essential point of the calculation because the coupling between p-p 
and p-h channels is appears mostly through interactions that have other than just leading-logarithmic flow. This is 
a special feature of the square or almost square Fermi surface and can be handled only by the non-local (outer-shell) 
contributions to the flow, using the Polchinski equation. 

The first aim of our calculation is to find the temperature at which the system flows towards strong coupling. 
We associate this temperature with a mean-field like critical temperature and call it T c flG . A typical mean-field 
theory then is regularized for T < T c by adding counterterms that contain fermions bilinearly coupled to some order 
parameter. In our theory the order parameter is not known a priory: it is determined by the manner in which 
the function U(9\, 82,83) diverges at T = T^ G . We perform a detailed analysis of the behavior of the angle-resolved 
correlation functions for antiferromagnetism and superconductivity and obtain the type and the symmetry of the order 
parameter determining the dominant correlations near the T^ G . The final result is the analogue of a mean-field phase 
diagram of the Hubbard model. We are considering a two-dimensional system where one should be careful about the 
interpretation of T c : in the case of magnetism, this indicates the onset of well- defined finite-range correlations. For 
weak interactions, this is typically a very well-defined crossover .0 In the case of pairing T^ G can be identified with 
the onset of quasi-long-range order. However, in real systems like copper oxides even a weak inter-plane two-particle 
hopping (particlc-hole-pair hopping for antiferromagnetism or Josephson tunneling for superconductivity) stabilizes 
a 3D long-range order. 

In section II we begin by the formulation of the many fermion system on a lattice in terms of functional integrals. 
We introduce the concept of the effective action and show how it can be formally calculated using the partial trace 
technique. We than derive the Kadanoff-Wilson-Polchinski exact renormalization group as one possible strategy for 
calculating the effective action in terms of the renormalization group flow of all vertices. Truncating the effective 
action at the level of sixth order vertices we obtain the one loop renormalization group equations for the effective 
interaction and for the selfenergy. In section III we apply the zero-temperature one-loop renormalization group to 
the Hubbard model on a square lattice. We derive the flow equations for the effective interaction function and for 
the angle-resolved correlation functions of superconducting and antiferromagnetic type. After discretization of the 
angle 8 on the Fermi surface we integrate numerically the flow and present the resulting phase diagram. In section IV 
we introduce finite temperature explicitly in the renormalization group equations. We then calculate the fixed point 
values of the correlation functions at temperatures near the instability. The conclusions are given in section V. 
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II. FORMULATION OF THE RENORMALIZATION GROUP FOR A MANY— FERMION PROBLEM ON 

A LATTICE 



The simplest model for interacting fermions on a two-dimensional square lattice is the Hubbard Hamiltonian 
H = -t ( a i* a 3,<? + a j,a a i,<r) + Y^2 UlUl "^^2^1 (2-1) 

where ai iCr (a| cr ) is the creation (annihilation) operator of an electron at the site i with spin a, t is the inter-site 
transfer integral, /i is the chemical potential, and Uq is the on-site Coulomb repulsion. After Fourier transform, the 
Hamiltonian writes 

H = Y1 &4k°<* + 2 Yj C/ " a -<T,k 1 +k 2 -k3 a -^ 2 alk3 a ^i » (2-2) 

crk er ki,k2,k3 

where 

£k = — 2t(cosk x + cos k y ) — /i (2.3) 

and the momenta are within the first Brillouin zone. In this section we want to derive the renormalization group for 
a more general problem. For that purpose we allow £k to have a general dependence on k. Furthermore, we suppose 
that the interaction can be nonlocal and dynamical, that is, we suppose that it depends on energies and momenta of 
the interacting particles. 

The statistical mechanics of souch general model is given by the partition functiorS 

Z = I V^!V^ e s{ * } (2.4) 



where the functional integration is over Grassmann variables 4 f ($) for all electrons in the Brillouin zone. The action 
S is given by 

+ \Y. T3 UoiK^K^K^K^K^K^^. (2.5) 

aa 1 u ni ,U) n2 ,u„ 3 ki,k2,k3 

The variables \E , (4') are labeled by the energy-momentum vector K = (cj„, k). £t is the bare spectrum measured from 
the Fermi level: 

£k = e k - fi , 

where e k is the band dispersion and [i the chemical potential. The energies and momenta are conserved so that 
Ki{Ki, K2, K3) — (u> ni + Lu n2 ~ w n 3 , ki + k2 — ^3). Uq{K\, K2, K3) is the most general spin independent interaction, 



a function of the frequencies and momenta. The derivation of the action (2.5) for a general model is equivalent to the 
derivation for the Hubbard modeBS, provided that we put Uq(K\, K2, A3) instead of the constant Uo and keep £ k 
general. 



G 



We want to derive the low energy effective action (LEEA) for this model. The low energy modes are the electronic 
degrees of freedom close to the Fermi surface. We will use this criterion and use the energy variable to discriminate 
fast (high-energy) modes from the slow (low-energy) ones \E'<. Let's choose some arbitrary nonzero high energy 
cutoff A defining a shell of wavevectors around the Fermi surface. The electronic variables can then be written 

V<r,K = 0(|£k| - A)¥> )(y , x + 9(A - |£k|)*<,a,K. (2.6) 

The slow modes are inside the shell ±A while the fast ones are outside, with |£k| going up to the physical cutoff Ao 
taken to be equal to the bandwidth so that we are sure that the whole Brillouin zone is taken into account. Note 
that the cutoff is imposed only on momentum space, while the Matsubara frequencies remain unlimited. The LEEA 



>S'a{ v I'<} is an action containing only slow modes and gives the same partition function as S eq.(2.5), or formally 



Z = / X>* < OT < e SA{ * <} , (2.7) 



This means that S'a{ v I'<} is calculated by taking the partial trace over only fast modes in eq.(2.4) 



£a{*<} = In J X>*>X>*> e s{ *<^> } (2.8) 

The LEEA contains a new effective kinetic part So a with a finite selfenergy term and a new interaction Si a- We have 
chosen that the Fermi surface for the bare electrons plays the role of the zero energy manifold. This still does not 
mean that we make the approximation of a Fermi surface unrenormalized by interactions: even if the Fermi surface 
of the LEEA is different from the bare one we are still allowed to use £g for the bookkeeping of our mode elimination. 



If we consider the slow modes as parameters, expression (2.8) can be evaluated, at least formally, using the linked 
cluster theorem.il The result is composed of three terms: 

s A {*<} = S'{* < } + n > + <y5{* < } (2.9) 

Only the interaction part Si of the action S can mix slow and fast modes: 

Si = + S/{* >} + Si{^ < , (2.10) 

while So is diagonal and can contain only one kind of modes in the same term: 

So = $>{*<} + So{*>} • (2.11) 

The first term in the equation ( [2.9] ) is then only a constant from the point of view of the fast electrons and is equal 
to S'o{^ r <} + S'j{^ , <}. £1> is the grand potential (times f3) of the fast electrons as if they were decoupled from the 
slow ones: 

fi> = -^ln(l + e~ /35k ) + ^(all connected clusters with Sj{*>}) . (2.12) 

k> 

This term gives only a shift of the total free energy of the system. 



The term <5S'{\E f < } in cq.(2.9) is the most interesting one. It brings the corrections due to the scattering processes 



of the slow modes on the fast ones into the LEEA and is given by the sum of all connected graphs composed of 
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FIG. 1. A few lowest order cumulants for <5S'{^'<}. All internal lines are integrated only over the fast (>) modes. 



S/{ , i , >} + S'/{ , I'<, V E'>}. If we draw the slow modes as external legs, the diagrams for (55'{\E , <} are the clusters with 
at least two legs. A few low order diagrams for 5S'{ V I'<} are given in fig.0. The terms with two external legs, labeled 
by a, b, c, and d in figure |l] are the selfenergy terms, renormalizing Sq. The terms with four legs: e, /, g and h 
renormalize the quartic interaction term Si. The terms with six (i and j) and more legs are new! They are created 
by the mode elimination procedure. 

Of physical interest is the LEE A for the electrons in the very vicinity of the Fermi surface (A <C Ep). Even if the 
coupling is small some of the loop diagrams will attain at low temperature (T < A) large values depending of the 
form of the Fermi surface. For example, if the Fermi surface is not close to van Hove singularities the particle-particle 
(p-p) diagram (/ in fig. |l|) for the four-point vertex with zero center of the mass momentum always has a logarithmic 
dependence like log(Ao/A), where Ao is the initial cutoff equal to the bandwidth. If the Fermi surface is nested the 
particle-hole (p-h) diagram (e in fig. |l|) at 2kp behaves in the same way. In the Hubbard model close to half-filling the 
van Hove singularities make both loops squares of logarithms. The perturbative calculation of the expansion ( |2.12 ) for 
small Uq is thus not straightforward: at least some of sets of diagrams, containing both p-p and p-h subdiagrams, have 
to be summed entirely. The lowest order diagram of that kind is the one denoted by h in fig. |} On the other hand 
the truncation of the LEEA at fourth order is in general allowed for weak coupling. However the direct summation 
of cumulants for J5{^<} (like T-matrix or RPA summation) can be performed in a useful and controlled way only 
for a limited number of physical problems, that is when some subsets of diagrams are dominant. The direct parquet 
summation for a general Fermi surface in more than one dimensions is probably very hard. It has been done only for 
the case of perfectly nested (fiat) Fermi surface@@. 

The problem is even more difficult if the coupling is not small. Then the criteria of most important sets of diagrams 
are not clear any more and even the truncation of LEEA at quartic or sextic term in ^< is not justified any more. 
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FIG. 2. The division of the Brillouin zone into the outer-shell (>), the on-shell (I), and the slow (<) modes. 

A. Kadanoff— Wilson— Polchinski renormalization group: Exact formulation 

A tractable way to construct the exact LEEA is to use the Kadanoff-Wilson-Polchinski renormalization group. 
Let us call the initial cutoff (the bandwidth) Ao and parameterize A by the renormalization parameter / so that 
A = AoCxp(— I). The idea of the renormalization group is to consider the transformation S = Sa - * ^A a exp(—l) as 
an infinite set of infinitesimal mode eliminations 

Sa — > ^A exp(-dO "~ * ^A ) cxp(-2d/) ~ * — — * ^A oxp(-i) 

At each step we eliminate Adl of modes at a distance A from both sides of the Fermi surface. We will see that the 
mode elimination of an infinitesimal shell of degrees of freedom is much simpler than the one-step procedure discussed 
in the previous section. 

From now on we will call the LEEA simply the effective action because in the process of successive mode elimination 



(2.13) A can have any value between Ao and zero. Indeed, it is of physical interest to follow the flow of the effective 
action S\ as A decreases. 

We now concentrate to one single step I — ► I + dl of the mode elimination. We call outer shell modes the modes 
already eliminated by the previous steps (the fast (>) modes) . The modes inside the shell [Ao exp(— I)— Aq exp(— I— dl)] 
are the ones to integrate out. We call them on-shell modes and denote them by (I). Figure || shows the division of 
the Brillouin zone into three types of modes (>, I, and <) for the case of the non-half-filled Hubbard model. 

If I is not the very first step, the effective action Sa contains couplings of all orders. Schematically it reads 

Sa = Sqa + Sia = r^ ** + ri°M*# + rfffww + ... (2.14) 
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The summation over all frequencies, momenta, and spins is assumed. The two point vertex T2 defines Wick's theorem 
at step I. In particular the propagator of the "bare electrons" T^ 1 changes as we proceed with the renormalization. 

The construction of the effective action one step further (at I + dl) is of the same form as the equation (|2.S|), 
with the difference that now the on-shell modes play the role of the fast modes. As we are interested only in the 
renormalization of vertices, we can skip the constant Q and we get the recursion relation: 

SA(i+di) = S A( i) + SS(l) . (2.15) 

The contribution 6S(l) is due to the elimination of I modes. It is given by the sum of all cumulants made of Sia{^i} 
and Sja{^ , !>^<} with two or more legs but now with all internal momenta constrained to be on-shell. Now we 
use the fact that dA is infinitesimal: In the expression for 5S(l) only the terms linear in dl will survive, to make 



the recursion (2.15) a differential equation for S\m. Generally the cumulants with m internal lines are proportional 
to dA m because every internal line is constraint to the shell. In principle only diagrams with one internal line are 
proportional to dl. If we group the terms with equal number of legs, we obtain the flow equation for vertices r„ , 
known as Polchinski equation for the vertices 

Si Only two types of diagrams with one internal line are possible: 
tree diagrams and loop diagrams. The Polchinski equation for the vertices is shown on figj^. Its symbolical form for 
a two-point vertex (2n — 2) is 

■£-T®(K) = -TJ2 [ d 2 fcTf(^',^,X',if)G ; (^')- (2-16) 
c A i ~r JdK 

This means that only the loop term renormalizes the selfenergy (see figure |^) . Both loop and tree terms are present 
in the Polchinski equation for the higher order vertices: 

a 



h,h JdA 

- TJ2 [ d 2 kT^ n+1) (K,K ll ...,K n ,K,K n+ll ...,K 2n )G l (K) . (2.17) 

The two-point vertex defines the one-particle propagator Gi at each step of the renormalization: 

G l (K) = (T«\K))- 1 (2.18) 

We use this renormalized propagator to construct the Wick theorem. We name T^^Ki, K n , K n +x, K2 n ) the 
vertex with 2n external legs at the step I of the renormalization, with legs {K\, .., K n } coming in and {K n+ i, i^2n} 



coming out. Symbols I\ and I2 in equation (2.17) are disjoint subsets {I\ PI 1% = 0) of the energy-momenta such that 
Ii U I2 = {K\, K2n}- The sum runs over all such sets. We have skipped spin indices for simplicity. 

We see that the Polchinski flow equation is a functional equation because all vertices are renormalized as functions 
of momenta and frequencies. It gives the exact renormalization group flow of the model. The careful reader will 
perhaps hesitate at this point: some of the loop diagrams with two internal lines as for example the p-p diagrams are 
also proportional to just dl if we put the total momentum and energy to zero. The same "anomaly" happens if the 
nesting is perfect for some p-h diagrams. In the usual renormalization group calculations only these contributions are 
taken into accountS because they give the dominant logarithmic part of the renormalization group flow. Then one 
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FIG. 3. Polchinski equation for the vertices with n = 2 and n > 2 legs. 



misses all non-logarithmic or "almost logarithmic" physics. This is not a problem in one dimension, for example. But 
if we are in two or more dimensions and especially if the nesting is good but not perfect, it is better to consider the 
exactly logarithmic configurations of energy- momentum as exceptions. If we formulate the flow equations correctly 
for the general case, the exact logarithmic terms will also appear, as we will see in the next subsection. 

In principle vertices of all orders are created with increasing powers of the initial coupling Uq: It is easy to see that 
the vertex I^n, (n > 2), is created by the tree term of the Polchinski equation with power (n— 1) of the bare coupling. 



This means that the truncation of the expansion (2.14) is equivalent to weak-coupling perturbation theory. 



B. Truncation of the Polchinski equation: one— loop renormalization group 

The one-loop renormalization for the vertex T4 (or for the effective interaction U{) is the perturbative procedure to 



truncate the flow equations at order U 2 . All terms of order higher than six in the expansion (2.14) are created with a 
power higher then 2 of the interaction by the tree term of the Polchinski equation. Thus putting Tg = Tio = ... = 0, 
we generate the one-loop renormalization group. The only contribution to the vertex Tg is than the tree term, made 
of two Ti terms connected by one line (see figure fja)). The line denoted with I has in principle to be taken dressed 
by the selfenergy at the step I defined as 

s ; = - r° = rf - iu n + £. (2.19) 

We assume that the selfenergy remains diagonal upon renormalization. This is consistent with the weak-coupling 
treatment because off-diagonal terms would imply the existence of some form of long- range order, which is out of 
the reach of the present calculation. All we can possible expect from our calculation is a divergence of some effective 
interaction signaling the onset of long-range order. 
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FIG. 4. (a) The six-point vertex for one-loop renormalization group, (b) Relation between the vertex T4 and the interaction 
U . (c) Recursion for the one-loop renormalization of the interaction U . 

We now go back to the formulation in terms of the interaction as denned by eq^^) and illustrated in fig.^(b). We 
will skip the spin indices where they are not necessary. The differential flow of the six point function at step I is 
according to the figj||(a) given by 

dT^(K 1 ,K 2 ,K 3 ,K 2 ,K^ = TY j [ d 2 k <5(k - ki - k 2 + k 3 )<5. 

. . JdA 



J U> n — UJ n l — UJ n 2 + U) n 3 



(2.20) 



The phase space integral is over the shell of thickness dA corresponding to step I, and Gi (K) is the renormalized Green 
function at the same step. Physically this is the propagator of an on-shell electron renormalized by the scattering on 
the fast electrons. Ui{K\, K 2 , -K3) is the effective interaction at step /. The vertex Tq at some step I is the integral of 
eq.( 2.20| ) over all steps between A = Ao and A/, that is over all fast degrees of freedom. On the other hand, there is 
no loop integration in this term: the Dirac function in equation ( [2.20 ) reduces the integral f d ^d 2 k to a single point 
k = ki + k 2 — k 3 (= k3/ + k,j' — k 2 '). To get T§ we can thus skip the integration over dl and take just care of 
momentum conservation. The effective action at the step / then reads 



Si 



^ ^K,cr(iU!n 



•& + Ej(JQ)**,ae(A(0-|&|) 



^ 3 E E ^(^i^2,^3)e^S )k3 , k4 ^ 3 ^x 4 ^x 2 ^ a 



aa 1 1,2,3 

+ T 5 E [<!: 

a.a' ,a" 1,2, 3, 2', 3' 



^ (A S, k3 , k ' k ' k '0(ia|-A(O) 



(2.21) 



where I' = In Ao/|£ k | (i. e. £k = A(Z')) is the scale fixed by external momenta, K = K\ + K 2 — K3 and the energy- 
momentum 4' = 1 + 2 + 2' — 3 — 3'. The summations over 1,2,3,... run over corresponding Matsubara frequencies 
and momenta. The term of the sixth order contains the interactions and Green functions from former steps I' < I of 
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the mode elimination since only fast degrees of freedom contribute to T l 6 so that I' < I. This constraint is imposed 
to the scxtic term of action ( 2.21 ) by 0(|£k| — A(Z)). The functions S^^k 2 k 3 k 4 an< ^ ®k^k 2 k 3 k' k' k' constrain the 
momenta in arguments to be slow modes (inside a shell of thickness ±A around the Fermi surface). This simply 
means that the fields described by the effective action at cutoff A(^) are inside the cutoff range. 

If the initial interaction Uq is spin independent, the renormalized interaction Ui will remain spin independent as 
well. It is thus not necessary to worry about spin indices and all two-particle interactions are given only by one 
function Ui(l, 2, 3). The detailed justification for that is given in the appendix The differential flow of £7j(l, 2, 3) is 
readily obtained applying the loop term of the Polchinski equation (the second term in equation (p7|) and in figure 
||(b)) to the six-leg part of the effective action ( |2.21 ). At first sight, two kinds of diagrams are created: one-particle 
reducible (1PR) and one-particle irreducible (1PI) ones. We will show that only 1PI diagrams contribute to the 
renormalization of the effective interaction : We can try to construct the 1PR diagram by contracting legs 2' and 4' in 
figure ^(a). This immediately implies that the internal line denoted with I and the leg 3' carry the same momentum. 
This momentum corresponds to some fast mode since line I is already integrated out. The conclusion is that the 



resulting four-point 1PR vertex can not be a vertex of the effective action ( 2.21 ) since this action contains only slow 

modes. Consequently, only 1PI diagrams renormalize the effective interaction between slow electrons. The resulting 

diagrams are all topologically different two-particle loops as shown on figure ^(c). The first diagram is a p-p diagram 

and the others are p-h diagrams. Let us ilustrate how we obtain the first diagram in fig. |](c). The procedure is shown 

on fig.|. The dia gram represents the p-p contribution to the effective interaction {7;(1, 2, 3) due to the elimination of 

the infinitesimal shell at step I. We take the six-leg diagram with the configuration of external momenta shown in the 

figure and with legs K being on-shell. Their contraction (dashed line) is done precisely at step I. The contraction K' 

was done at a previous step l pp fixed by momenta K, K\, and K<i (see equation ( 2.20D ). K\ + K2 = Q pp = {ui npp , q pp ) 

is the total energy-momentum in the p-p process. The scale l pp is then given by 
t 

l PP = - ] n^r^- (2-22) 

Similar constructions give all other (p-h) diagrams. One has to take care of both direct and exchange interactions 
(fig. ^(b)) to get four different graphs. The interactions and one-particle propagators are to be taken at the scale 
l p h < I. This scale is determined by the momentum transfer q p h as 

lph = _ ln &+^H (2 23) 

In second, third, and fourth diagram in fig. |](c) the energy-momentum transfer is Q p h = K\ — K3, while in the last 
diagram Q p h = K± — K4. 

We see that even though the Polchinski equation appears as local in I, the flow at step I depends on the Green 
functions and interactions at the former steps l pp ,l p h < I- The reason for that is the dependence of the six point 
function on two-particle interaction and on one-particle propagator at all former steps. 

In fig.[|(c) the internal lines labeled with I are to be integrated over two infinitesimal shells of the width dA at 
£0 = iAj. For this purpose we pass from the Cartesian measure dk x dk y to the measure 



13 



K 



A K' 



K 



FIG. 5. Construction of the p-p diagram from the six-leg vertex. 



where s are lines (or surfaces in D > 2) of constant energy £(k) and v(s,£') is the group velocity as defined by 
9£(k) / dk± (k± is the component of the momentum perpendicular to the equal-energy lines. We will use measure 
( |2.24 ) in what follows, where we write the analytic expression for the flow of function E7j(l, 2, 3). 
From the diagrams in fig.^(c) one obtains the following expression 



ai 



(3 PP {U, U} + 2(3 ph {U, U} - P ph {U, XU} - P ph {XU, U} - Xp ph {XU, XU}. 



(2.25) 



/3 is a four point object and a bilinear functional of Ui>, {V < I). The operator X is the exchange operator acting on 
a four-point function: XF(1, 2, 3, 4) = F(2, 1, 3, 4). j3 pp and (3 p h are the p-p and p-h parts of the (3 function given by 



f3 pp {U,U} = (~{U,U} + E{XU,XU}) 
Pph{Ui,U 2 } = {Il{Ux,U 2 } + TIi{U u U 2 }) 



(2.26) 
(2.27) 



T is the time inversion operator acting on a four point function: TF(1, 2, 3, 4) = F(3, 4, 1, 2). The functions 3 and II 
correspond to the on-shell integrals of the p-p and p-h bubbles: 



= E /^e(l&- te l-Ai)rX;G«(^ M )G w (-^M + O w ) 



( 2 -) 2 „ 



(2.28) 



^{U^U 2 }{K U K 2 ,K^K^ 
-A, 



—© (ia„ +qp J - A ; ) Gi{K (u) ) G lph {K {v) - Q ph ) . 

u=-r,— uj n 



x U Uth {Kx,K {v) ,Ks)U*,i th (K 4 , K {v) ,K 2 ) 



(2.29) 
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Ui et U2 can be U or XU as required by eq.(2.25). The summation over index v — +, — is over two shells at £0 = ±A;; 
the velocities are v v — v(s u ,t; — isA) and K v symbolizes (k„,w„). The quantity 



is the total energy-momentum and 



Qph = {u n ,ph,<lph) — Ki — K 3 , 



is the energy-momentum transfer between the currents (1,3) and (2,4) where 1,2,3 and 4 are the external variables of 
S and II. The scales l pp and l p h are defined by expressions ( p. 22 ) and ( p. 23 ). As already discussed, l pp and l p h depend 
on the integration variable k„ and on the configuration of the external energies-momenta. Let's call the external legs 
of the total (3 function ( 2.25|) 1,2,3 and 4. Note that the operator X exchanges the external legs 1 and 2 in the last 
term of this expression. This means that the energy-momentum transfer in this term is Q p h — Kf — Kg and not 
Ki — as in the first three e-h terms. In the standard language (see for example ref. ^) the p-h terms with transfer 
1 — 3 are called zero-sound (ZS) terms, and the terms with transfer 1 — 4 are ZS'. 

Let us explain briefly how we obtained the flow equation ( |2.25 ). The first term is simply the p-p loop with 
U interaction. The remaining terms are different versions of the p-h loops, corresponding respectively to the p-h 
diagrams in the figure |](c). They can all be seen as a single loop (3 p h (with the topology of the second diagram in 
figj^(c)), given by (2.27) and ( p. 29 ) by performing appropriately the operation X. The third and the fourth graph 
can be drawn as the second one by one exchange: In the third term we replace the upper interaction line U by XU 
and in the fourth the lower one. After this manipulation both diagrams look like the second diagram. The last graph 
is more complicated: one has to perform X upon both interactions and upon the whole graph to see it as f3 p h- The 
factor 2 before the second term is due to spin summation in the loop. All other diagrams have fixed spin. 



The flow of the effective action ( 2.21 ) is still not completely determined because we do not know how the self-energy 
Y,i(K) is renormalized. The differential flow for ~Ei(K) is readily found from the Polchinski equation for the two-point 
function (2.16) shown graphically in figjU In the language of the effective interaction £/;(l,2,3) this gives Hartree 
and Fock like contributions shown in the figJJ(c). We get the renormalization equation 

J^TT = aHartree{f//}(ifl) + OFbck{^}(^l) + "JTMO- (2-30) 

The first term is the Hartree term 



a H artrec{t//}(^l) = 7^ E / G/ (iC„) i(3 - X)Ul(Kx, K v , K{) (2.31) 

and the second is the exchange term 

a F ock{C//}(^i) = E f^-T^G^K^fl-XpiiK^K^K!). (2.32) 

The third term is added to cancel the chemical potential renormalization due to the homogeneous part of the direct 
term 

<° m (0 - --^2 / rf2fc l I>Hartree{^X*Tl) • (2-33) 
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If the initial interaction has no dynamics, the first nontrivial contributions to the flow of the selfenergy comes from 
the renormalized and not the bare interaction. As the interaction is renormalized by one-loop processes this implies 
that the interesting part of the flow of T,i(K) is given by two loops. On the level of the present one-loop calculation 
it is thus consistent to neglect selfenergy corrections. This is what we do in the subsequent one-loop renormalization 
of the Hubbard model. 



III. RENORMALIZATION GROUP FOR THE HUBBARD MODEL 

In this section we will apply the above renormalization-group procedure to the Hubbard model. The model is given 



by equations ( |2.l|) , (2.2), and (2^3). The initial action S for the Hubbard model is given by expression (|2.5| ), with the 
dispersion ( |2.3| ) and with the initial interaction Uq = constant [K\, K 2 , K3). The interaction will depend more and 
more on {K\, K 2 , K3) as we go on with the renormalization (as I increases) so that we will see at work the functional 
aspect of our renormalization group. We will complete our analysis by the renormalization of two-particle correlation 
functions. 



A. Renormalization of the interaction 



If we neglect the selfenergy corrections, the flow of the effective interaction is completely determined by the expres- 
sions (2.25- [2~29| ) with the bare propagators instead of the renormalized ones: 



Gi(K) -> G (K) = (iLu n - (3.1) 
The effective interaction is a function of three energy- momenta. This makes formulae ( |2.28 ) and ( 2.29|) very compli- 



cated. For that reason we will consider only the marginal part of the dependence of Ui on energy- momenta. This 
approximation is justified by the zero-order scaling and power-counting argumentsB@@ For example, in one di- 
mension this procedure justifies the well-known g-ology modelEl: one adds an index i to the electrons so that all 
electrons moving to the left have i = — and all right-movers have i = +. Then the marginal interactions do not 
depend on impuls k and energy u> of the electrons in interaction, but only on their indices i. This can be seen as 
parametrization of the interactions as if the electrons were on the Fermi surface (or points in ID) with uj = 0. In two 
dimensions the marginal interactions depend only on polar coordinates of the wave vectors. Only the interactions 
between electrons at the Fermi surface are then kept and, if the Fermi surface is not nested, one gets the LEEA for 
the Fermi liquid.SE3. The marginal processes in that case are 

V(9 1 ,9 2 ) = U 1 (K 1 ,-K 1 ,K 2 );cj 1 , 2 =0; £k lt2 =0 (3.2) 

and 

F(6 1 ,9 2 ) = Ui{K 1 ,K 2 ,K 1 );ui, 2 = 0; £ kl , 2 = 0. (3.3) 

V is the pairing amplitude and F is the forward scattering related to the Fermi liquid parameter. Both kinds of 
processes are shown in figure || for the case of the Hubbard model far from the half filling. We have analyzed in detail 



1G 



31 41 




V (67,63) 




F(e 7 ,e 2 ) 



FIG. 6. The marginal interactions in the BCS regime. 
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FIG. 7. Some of the marginal processes for the square Fermi surface. 

this problem in a former article. @ Let us now concentrate on the square Fermi surface, for the half filled Hubbard 
model. The processes between electrons on the Fermi surface are now labeled with three variables instead of two 
as in the case of the Fermi surface without nesting : if we put particles 1, 2 and 3 anywhere on one side or on two 
opposite sides of the square, the fourth falls exactly on the square as well. This is due to the perfect flatness of the 
Fermi surface and to the marginality of the umklapp processes. A few examples of marginal interactions between the 
electrons on the square are shown on figure ^. The interaction depends only on the positions of the particles on the 
square. The "angle" 9 can be defined in a way shown on figure |[ It is important to notice that even if the filling 
is not exactly one-half (and the Fermi surface not exactly square), all above interactions will still be important, as 
long as the effective phase space is open i.e. when A > (as on figure |J). We thus take as marginal all effective 
interactions viewed as functions of three angles 8 of the particles: 



Ui(K 1 ,K 2 ,K 3 ) -> Ui(6i, 9 2 ,6 3 ) ;wi,2,3,4 = ; 1,2,3,4 are on the square. 



(3.4) 
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FIG. 8. The organization of the mode elimination. Dashed lines are the lines of constant "angle" 9. 

When the cutoff becomes smaller than the chemical potential, we are back to the non nested case, in which the 
functions V and F are the marginal interactions. They read 

Vi(6x,6 2 ) = 17,(01,01 +ir,9 2 ) , F l (6 1 ,9 2 ) = Ui(6 1 ,9 2 ,e 1 ) . (3.5) 



Altogether, for the half and almost half filled Hubbard model the function C7, (6*i , 6*2, 6*3) given by (3.4) contains all 
marginal scattering processes. The renormalization group analysis is now much simpler because we deal with a 
function of three variables instead of six. 

We will now derive the flow equation for Ui{6\, 8 2 , 9 3) at zero temperature. The pleasant aspect of the Kadanoff- 
Wilson-Polchinski mode elimination technique at T — is that A can then be interpreted as the temperature. Namely, 
the interaction at some temperature T is renormalized mainly by virtual processes involving "quantum" electrons, 
those with energy larger than T, having almost the distribution of the T = electrons. This is exactly what we do 
with the renormalization group: only modes with |£k| > A are involved in the virtual processes renormalizing Ui- 
Consequently, A is not only the measure of how many electrons are already integrated out : it has a physical meaning 
of the effective temperature. 



Replacing Gi(K) by Go(K) and Ui{K\, K 2 , K 3 ) by Ui{9\, 9 2 , 83) in expressions ( 2.28 ) and ( 2.2S ), we can get inter- 
actions out of Matsubara summations and perform the summations analytically. After taking the T — > limit we 
obtain 

E{U,U}(9 1 ,6 2 ,9 3 ) = -^ J2 I *) 6 ^-^ 9 " A) x 

x [7^(0!, 2 , 0)£/ ipi ,(03, 04,0) , (3.6) 

n{ f / 1 ,t/ 2 } ( 1 ,0 2 ,0 3 ) = 7 ^ e /^A,0) e( '^;^^ (l ^ +q - J ' A) x 
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FIG. 9. The parquet summation. 



x Uu ph (e 1 ,9,9 3 )u 24ph (e i ,0,0 2 ) 



(3.7) 



k„ is the momentum of a particle at the angle with energy £ = vh. J{e,0) = J[(x,y)/(e,0)] — ds/d0/v(0, e) is the 
Jacobian of the transformation from rectangular coordinates in momentum space to polar coordinates. One should 
not forget that the scales l vp and l p h, which make the flow equation non-local, depend on external momenta and on the 
integration variable through q pp , q p h and k„ as given by the relations ( 2.22 . [2.23| ). U\ and U 2 represent U or XU as 



required by equation (2.25). Equation ( |3.6| ) gives the leading logarithmic flow in the p-p channel for the configuration 
of momenta with q pp — while equation ([^7]) gives the leading logarihmic flow in the p-h channel only exactly at half 
filling for q p /j = (7r,7r). In the standard renormalization-group procedure@, only these configurations are taken into 
account. We see that in our formalism they are taken into account on equal footing with all other scattering processes, 
with any values of q pp and q p ^ , as illustratd on figure 0. The processes with the leading logarithmic renormalization 
in one channel and with less strong but still important flow in the other channel (as the processes in figures 0(d) 
and (e)) are the processes which couple strongly both renormalization channels. For exemple, the process in figure 
0(d) has leading logarithmic renormalization in p-p channel and less strong (but still logarithmic!, because of partial 
nesting) renormalization in p-h channel, while the process in figure 0(c) has perfect nesting and, consequently, leading 
lorarithm in p-h channel and weaker logarithmic flow in p-p channel. 

If we want to see which series of diagrams is generated by our renormalization group we have to solve the differential 



equation ( 2.25 ) for Ui in iterations of the bare interaction Uq. The obtained series is exactly the parquet summation. 
It is constructed from all iterations of five basic loop diagrams from figure |](c) . A few lowest order parquet diagrams 
are shown on figure |^. 

An important aspect of the non half filled Hubbard model is that it can not be solved by a scale invariant renormal- 
ization group. The finite chemical potential determines the intrinsic scale, the physical interpretation of which is the 
crossover between two different renormalization regimes. The crossover can be seen from the explicit scale dependence 
of p-p and p-h differential loops and II(Z). We can define the quantities 



Pph(l) = n /{l, l}q ph =(7T^) 



(3.8) 
(3.9) 



19 



2.00 




FIG. 10. The quantities f3p P (l) and f3p h (l). The crossover is at I — Z M = 6 

They measure respectively the dominant parts of p-p and p-h renormalization tendencies. The configurations of 
momenta are chosen to give the most important flow : for the p-p channel at zero total momentum and for the 
p-h channel at the antiferromagnetic wave vector. The quantities 0p p (l) and (3p h (l) are shown in the figure [l^ for 
finite chemical potential /i = — Aoea;p(— /^). We see that for I < both differential loops have linear dependence in 
logarithmic variable 1; the total (integrated) loops are thus square logarithmic, as is known for the half filled band. 
When I > I ft the function (3p p (l) crosses over to constant which gives the logarithm of the Cooper bubble. f3p h (l) decays 
exponentially as exp(— 21) ~ A 2 : the nesting does not exist any more and the p-h flow crosses over to irrelevance. 
We call the first regime the parquet regime because both loops are important. The second regime, in which only the 
Cooper channel flows, we call BCS regime. The topology of the effective phase space in the parquet regime is open 
(see figure ^) and, in BCS regime, the phase space is a regular closed ring around the Fermi surface, as on figure 
||. The peak of f3p p (l) at I = 1^ is the enhancement due to van Hove singularity. The peak does not exist in Pp h {l) 



because of the 0-function constraint in (3.7). As we will see later, the renormalization in the parquet regime will 
give rise to precursors of a strong coupling fixed point with dominant antiferromagnetic correlations while in the BCS 
regime only a Cooper-like instability is possible. 

It is difficult to read from the sole flow of the interaction Ui , what kind of correlations are enhanced and possibly 
divergent. For that purpose we have to calculate the renormalization of the correlation functions. 
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FIG. 11. The angle-dependent correlation functions. 

B. Renormalization of the correlation functions 

It is well knownElIlllS that, in studying the anisotropic superconductivity, one has to consider the pairing 
amplitude as a function of two angles, V^Oi^O^)- The angles determine the angular positions of the Cooper pairs 
annihilated (#1) and created (# 2 ) in the scattering. This interaction will be intimately related to the superconducting 
correlation function \ (@i> $2)- I n the same spirit we can define the correlation function for the antiferromagnetism 
dependent on two angles. We will define both correlation functions in the following way 

Xq(0i,0 2 ;|n-T 2 |) 

= J de 1 J de 2 J(ei,0i)J'(e2,02)(A*(ei,0i;Ti)Aq(e 2 ,0 2 ;T 2 )), (3.10) 

with 5 = SC, AF ("superconductivity" or "antiferromagnetism"). The symbols ">" mean that the energy integrations 
run over energies outside of the shell ±A. Consequently, \ SC an d X AF are interpreted as the susceptibilities at the 
temperature T = A. They measure the response of outer shell electrons for given A. The order parameter variables 
are 

A* c (e, 9;t) = J2 ^,k(r)*- ff ,- k+q (r), (3.11) 

a 

A AF (e, 6- t) ee ]T * CT ,k(r)*_ CT , k+(7r , 7r)+q (r), (3.12) 

a 

where k is given by the angle 8 and the energy e. The figure |ll| illustrates what configurations of four angles are 
described by the correlation functions X SC {^1^2) and \ AF (^1,^2): the first measures the correlation of one cooper 
pair at 0\ with the other at 62 and the second represents the correlation of the momentum (tt, tt) p-h pair at 9\ with 
the other p-h pair at 62- The correlation functions can be seen as response functions of the system to an infinitesimal 
external field, as was done by Bourbonnais and CaronEl in one dimension. We will generalize this procedure to two 
dimensions. We add to the action Si = o the term 



S{h}i =0 = JdrJdqJdO J deJ(e,0)A q (e,0;T) 



h s {6; t) +h.c. (3.13) 
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FIG. 12. Tree diagrams conataining the source fields for AF. 



The external angle dependent fields h s (9; t) are the source fields coupled to the order parameter variables of the type 



8. The correlation functions ( 3. 10| ) are obtained as 

S 2 In Z 



6h s J9r,n)5M(62;T 2 ) 



(3.14) 

Putting slow modes to zero means symbolically that we want the response only from the fast modes, as defined in 



equation ( 3.10 ). We consider only the static and long-wavelength limit. For that reason we will simply write x(0i, 02) 
instead of Xq=o (01 , 02 ; = 0) and h(9) instead of h <i= o{9]iuj — 0). The correlations with the nonzero q and to are 
related to the dynamics of the collective modes, a problem which we do not study in this work. 



We now apply the Kadanoff-Wilson-Polchinski formalism to the action containing terms (3.13). The procedure of 



collecting differential cumulants is analogous to what we explained in previous section, but now we treat S{ h} terms 
together with the interaction part Si. To obtain the correlation functions for h — > 0, it is sufficient to follow the 
renormalization of the first two terms in powers of h in the /i-dependent part of the effective action. They read: 

r A(0 



S{h}i 



d6-L <p d9 



deJ(e,9 1 )A 6 (l (e,e 1 ;T) 



zf(9 1 ,9 2 )h d (9 2 ) +h.c. 



+ <* d$x & d9 2 h 5 (9 1 )xf(9i,9 2 )h s (9 2 ) + tree terms{/i/i} 



(3.15) 

The term with xf(9i,9 2 ) contains no electronic variable: it results from the elimination of all outer-shell electrons. 
From the definition ( |3.14[ ) one can see that xf(9i,9 2 ) is just the susceptibility of type S. The "tree terms" are the 
terms containing one outer-shell contraction, two slow-electron fields and fields h and h. They are illustrated by 
figure |l^ for the AF channel. The square symbolizes the effective interaction for antiferromagnetism V t AF . We obtain 
it from the spin-spin interaction U a (see appendix |X|) putting the particles 1 and 3 on the opposite sides of the square 
Fermi surface so that ki k3 = (±7r, ±7r) (as on figure |ll"| ): 



V l AF (9 1 ,9 2 )=-(XU)(9 1 ,9 2 ,9 1 ) 



9 is a function of 9 such that 



k(0) - k(0) = (7T,7T) 



(3.16) 



(3.17) 



k being on the square Fermi surface. The tree terms for SC channel are analogous, but with different orientations of 
the arrows: in the vertex zf c both arrows point outwards and in the vertex zl sc both arrows point inwards. The 
corresponding interaction is the familiar effective Cooper amplitude V ( |3.5| ). All "tree terms" in equation ( |3 .15 ) are 
produced by the tree term of the Polchinski equation applied to the action with the S{h}i terms. 
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FIG. 13. The recursion relations for the vertices and for the correlation functions of the superconducting and of the antifer- 
romagnetic type. 



The coefficient zf(6x, 9 2 ) is the effective vertex of type 6. The equation ( |3.13 ) gives the initial conditions for z 
zf= o (9 1 ,0 2 ) = 5 D (0 1 -6 2 ), (3.18) 



where 8d is the Dirac function and for x 
X? =o (0i,0 2 ) = O. 



(3.19) 



The differential flow of the triangular vertices zi and of the correlation functions x is obtained from the loop 
diagram of the Polchinski equation applied to the tree terms in the action S{h}i. For the AF channel the, cumulant 
with on-shell integration of the electrons A and B on figure |lj gives the contributions to the vertex zf F and the 
cumulant with electrons A' and B' on the shell contributes to the susceptibility Xi ■ A similar construction yields 
the renormalization of the vertex and of the susceptibility for the superconductivity. The resulting diagrams for the 
differential recursion relations for both channels are shown on figure The corresponding flow equations write 



if (#1,02) 



de z f s (e 1 ,e)D s l (0)v l s s (0,0 2 ) 



(3.20) 



and 



(3.21) 
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The scales l$c and Iaf symbolize the scales l pp and l p h given by expressions ( 2.22| ) and (2.23), with the total 
momentum q pp = and with the momentum transfer q p h = {-k, it) (antiferromagnetic wavevector): 



h = 



l PP \ci PP =Q = for 5 = SC 

Iphl^^ri = In a^2H for 5 = AF 

We see that the renormalization of the antiferromagnetic correlation function is non-local in I if the filling is not 
exactly one-half. The function Df(8) is 

Df C {0) = \ £ J(vA(l),8) (3.22) 
for the superconducting channel and 

A W -2 1 + | M |/A(0 ' (3 ' 23) 
for the antiferromagnetism where only the negative shell (y = —1) contributes to the flow. One sees that Df F (9) 
decays exponentially with I for A <C \fi\: in the BCS regime the correlation function for antiferromagnetism saturates 
with increasing I. 

From the equations ( |3.20 ) and (3.21) we see that information about the symmetry of the correlations is determined 
from the symmetry of the effective interactions : functions Df(0) have a total lattice symmetry, but the interactions 
Vf{6,62) can belong to any of the representations of the crystal symmetry group, in our case the D4 point group. 
The decomposition of the interaction in terms of all basis functions of all irreducible representations of the D4 group 
is discussed in detail in our previous paper.il The diagonalization of the correlation functions xf (#1,6*2) gives the 
final answer about which correlations are dominant in both AF and SC channels. The strength of the dominant 
correlations is associated to the maximal eigenvalues and the corresponding eigenvectors determine the symmetry and 
the form of the microscopic fluctuating field. 



C. Discretization of renormalization group equations 

The interaction Ui that we want to renormalize is a function of three continuous angular variables 9\, 62 and 



The beta function given by the equations (2.25-2.27, 3.6, 3.7) is a complicated function bilinear in U, and it does 



not seem possible to find an analytic solution for the flow of the interaction. We thus use numerical method. For 
that purpose we cut the Brillouin zone in m, angular (9) patches (see figure ||) and we assume that the interaction 
is a function only of the patch indices (ii, 221*3) of the three angles 0\, 62 and 63. After the discretization of the 
interaction function the differential loops 5 and II become also functions of three indices: 

rrii 

E{U,U}(ii, i 2 , 13) = ^2 Bpp(ii,i 2) i; l)Ui pp (i 1 ,i 2 , i)Ui pp (i 3 ,i4,, i), (3.24) 

j=0 

mi 

n{i/i, U 2 }(ii,i 2 , «a) = £ B ph (ii,i z , i; l)U Uph i3)U 2 ,i ph («4, i, i%) (3.25) 

i=0 
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with 

^(^0 = 7^ E / dBJ(uA,e f {v ^-^? f { ^-^- A) (3.26) 

and 

B vh (i x ,i z ,v,l) = > / ddj(vk,6,)^ '— ± '-. (3.27) 

The total momentum and the momentum transfer become discrete variables: 

q pp = k(ii) + k(i 2 ), 



q p h = k(ii) - k(i 3 ). 

The integral J*^, is over i-th angular sector. 

For a given number m, of patches the number of coupling constants is equal to the number of configurations of 
three indices for all four particles lying on the square Fermi surface. That is a very large number. However because 
of the symmetry many of the coupling constants are identical. The available symmetries are: (i) The symmetries 
of the Z?4 point group (mirror, 7r/4-rotations); (ii) Time inversion symmetry T, exchanging particles with holes and 
vice versa (see appendix ^); (iii) The exchange symmetry: it is allowed to exchange simultaneously (1,2) and (3,4) 
particles; (iv) The freedom of choice of the points at the edges of the Brillouin zone. Figure [l4| illustrates some of 
the symmetry operations applied to one of the coupling constants. The same figure shows the relation between the 
number of patches and the corresponding number of different marginal coupling constants, and the list of the coupling 
constants for m, = 4 and Wj = 8 . 

The renormalization of the interaction as a function of three angles is now represented by a set of coupled differential 
equations, one for each coupling constant. In the same way we discretize the correlation functions xf ^2) and the 



vertices zf (61,62). The equations (3.20) and (3.21) become 

zf(i 1 J 2 ) = -Y^4Ah,i)Df(i)V l s s (i,i 2 ) (3.28) 

i 

and 

*? (ixM) = - E 4 s {h,i)Di(i)z 5 h (i, h) (3.29) 

i 

with 

Dfii)^ [ d6Df(6). (3.30) 
J[i\ 

The initial conditions are the same as in the continuous case, provided we replaced the delta function by the Kronecker 
symbol divided by m,: 

5 D (6 - 6') ^ 5 iti , frm. 
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FIG. 14. The figure shows how we reduce the number of coupling constants by applying symmetry transformations. The 
dependence of the number of independent coupling constants on the number of angular patches m; and the list of coupling 
constants for m 4 = 4 and m, = 8 are also shown. 
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FIG. 15. The flow of a few typical (among 93) scattering amplitudes for a Fermi surface covered by 16 patches, for chemical 
potential \fj,\ = 8iexp( — 7.8) and initial interaction U = 4t/3. 

D. Results and discussion 

We have integrated numerically the renormalization equations for all coupling constants and for correlation functions 
and we have analyzed how the results change as functions of the initial interaction U$ and of the chemical potential 

We first look at the renormalization flow of the coupling constants. Figure [15] shows the flow of several (among 
93) coupling constants for a Mi =16— patches discretization; the choice of the input parameters is Uq — 4i/3 and 

= ln8i/|/i| = 7.8. The divergence happens at the critical scale l c w 5.3. Approaching this point, some of the 
coupling constants increase and diverge, while some decrease and, after changing their sign, diverge to — oo. Some do 
not change significantly upon renormalization and do not diverge. For example the coupling constant t/(0,mj/2, 0) 
diverges very strongly to — oo. It is a typical interaction with singular Cooper channel (q PP = 0) without nesting. 
Indeed, all coupling constants obeying only Cooper condition (q pp = 0) and without logarithmic flow in p-h channel 
diverge to — oo. This is what we expected since the p-p channel "pushes" interactions downwards in a repulsive model. 
However, instead of just decaying to zero, they continue to decrease towards — oo because the Cooper amplitude obtains 
attractive components from p-h diagrams in, for example, the D-wave channel. The coupling constants with nesting 
between particles 1 and 3 or 1 and 4 diverge to +oo. Among the interactions with nesting there are also umklapp 
processes like (0,0,4) or (2,2,10) in figure The processes without divergence are those without any logarithmic 
instability neither from the nesting nor from the Cooper logarithm. 

The critical scale l c depends on the initial interaction and on the chemical potential. We associate the cutoff 
A = AoCxp(— l c ) to the critical temperature T c flG . Figure |l6| shows T c flG as a function of the chemical potential 
calculated for m, = 32 patches (497 coupling constants). decreases rapidly but never really falls to zero : it 

becomes exponentially small far from half filling, the regime analyzed in rcf. p2[ Our numerical calculations show that 
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FIG. 16. The phase diagram. The solid line is the critical temperature T c and the dashed line is the temperature Tj 1 



this form is universal if one measures \i in units of the critical temperature at half filling Xj?: 



RG 



= T?xf 



(3.31) 



where / is the universal function visible on figure [l^. Thus T c depends on the interaction only through T ( 
8iexp(— Z°), where 
At 



(I 



0\2 



c 



U 



(3.32) 

The dashed line in figure 16 represents the critical temperature T^' IF that one 



C being a numerical constant, C ~ I 
obtains when taking into account only the last term of equation ( 2.25 ): X (3{XU, XU}. This is the "renormalization 
group" version of the RPA summation, equivalent to the mean field for the antiferromagnetism. We see now the 
main difference between the critical temperature in the mean field approximation and the result obtained with the 
renormalization group: in the case of weak doping, because of the destructive interference between channels p-p and 
p-h, T^ G is slightly reduced with respec to T C MF . The ratio between the critical scales l^ G and l^f F (associated to 
RG and MF critical temperatures) do not depend on the interaction. Its value at half filling is l ] f F jl^' G = 0.985, 
which is not far from the value 0.981, calculated by Dzyaloshinskii and Yakovenko using parquet equations.0 T^ IF 
disappears completely at some threshold doping. This means that the physical mechanisms which reduce T^ G near 
half filling, at higher doping enhance T^ G keeping it always non-zero. 

The straight line T — fi is roughly the crossover between the parquet and the BCS regimes. If the instability occurs 
in the parquet regime, both p-p and p-h correlations are strongly enhanced near T^ G . On the other hand, in the 
BCS regime only the p-p correlations are critical. To know wich which fluctuations are the most important at the 
instability, we need the renormalization of the correlation functions. However, there is a formal problem related to 
the fact that we are performing the renormalization at T = and associating the cutoff to the temperature: the 



renormalization equations for the antiferromagnetic vertex (|3.28j) and the correlation function (p.29|) (S = AF) are 
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FIG. 17. The flow of V c sc (dashed line) and of V C AF (solid line) for \fx\/(4t) =0(A); 0.00067(B); 0.0018(C); 0.0049(D); 
0.0081(E). 




retarded in I of a quantity I — In Aj ^2|^ . If the interaction V AF diverges at I = l c , the divergence of the function \ AF 
will be retarded. Since we can not go further than I = l c in the renormalization this divergence can not be seen in the 
present formalism. The cure is to work at finite temperature. In that case one performs the full renormalization, up 
to I = oo and the final fixed point correlation functions are the ones at the given temperature. This is the procedure 
that we use in the following chapter. 

However, in the zero temperature formalism on can still get some idea of what happens with different correlations 
at I = l c : The finite cutoff divergence of the correlation functions for SC and AF are determined exclusively by the 
divergence of the effective interactions V^ and V AF '. Furthermore the symmetry of the correlation functions is 
also brought only by the symmetry of the effective interactions. It is thus reasonable to assume that the dominant 
eigenvalue of the correlation function is driven mostly by the dominant attractive (negative) eigenvalue of the cor- 
responding effective interaction. From the renormalization of the interaction (the set of coupling constants) we can 



deduce the flow of the effective interactions V t sc and V AF as given by equations (BJ5) and ( 3.16 ). The diagonalization 
is straightforward because both interactions are matrices whose rows and columns are labeled by disretized angular 
variables. Let's call the most attractive eigenvalues of V^ and V AF respectively Vf c and V AF . Figure [l?] shows 
the flow of V c sc and of V C AF near the critical point as a function of ln[(A — X" c (/i)) / At) for several values of the 
chemical potential. The critical temperature T c (//) is adjusted for every value of /i. Solid lines represent the anti- 
ferromagnetic interaction V AF . The corresponding eigenvector belongs to A±. It is a standard S-wave. The dashed 
lines represent the flow of Vf c . Its eigenvector belongs to the B\ representation (D x 2_ y 2-wecve). Both coupling con- 
stants are always enhanced by the renormalization, which means that the correlation functions are always enhanced 
respectively to their value at Uq — 0. The possibility of the charge density wave instability is excluded: we have 
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checked that all eigenvalues of the charge interactions U c (see appendix |A]) at 2kp decay upon renormalization. The 
competition between the divergences of Vf c and V AF is clearly visible in the figure. At half filling the coupling V AF 
diverges faster than Vf c . As the chemical potential increases, both divergences are weaker but in V AF an inversion 
of the slope is visible. This is the signature of the crossover from the parquet to the BCS regime. At half filling this 
crossover does not exist and the slope of V^, is always upwards. The lines labeled by (C) in the figure correspond 
to the critical temperature in the parquet regime. However, V AF starts to "feel" the proximity of the crossover: the 
critical scale is l c = 6.025 and the crossover occurs at about l x ~ 1^ = 7. The divergence of V AF can still entail the 
divergence of the antiferromagnetic correlation function because the nesting is still relevant. The Cooper amplitude 
Vf c always has an upward slope and diverges at T c because the p-p channel has a logarithmic instability for any 
doping. Lines (D) and (E) are examples the flow in the BCS regime. After some saturation tendencies, V AF still 
diverges at Tf G . This divergence is only due to the p-p loop: for a choice of angles 9\ and 02 such that 62 = 9i + tt, 



the function V AF (0i, 62) = V sc {6\, 0i), that diverge. The relation between 9\ and 9\ is given by equation (3.17). The 
interaction V AF is thus driven upwards by the Cooper channel. It has no effect on the correlation function for the 
antifcrromagnetism because its flow has disappeared together with the nesting. This will become visible in the next 
section where we calculate the temperature dependence of both correlation functions near the critical temperature. 



IV. FINITE TEMPERATURE RG 

In the zero temperature formalism the flow of different quantities was of physical interest. In the finite temperature 
renormalization group, we are interested in the fixed point value of the correlation functions. The temperature is 
taken as input parameter. If T is larger than the critical cutoff A c (called T c flG in the previous section), the divergence 
of the renormalization flow will disappear. Consequently we will be able to control the flow all the way down to the 
fixed point A = 0. In the zero temperature formalism the effects of the elimination of the slow modes are neglected. At 
finite temperature all modes are integrated so that contributions of the thermal electrons are also taken into account. 
The other advantage of the finite temperature renormalization group is that we can explicitly follow the temperature 
dependence of the correlation functions for the superconductivity and the antifcrromagnetism. 

Formally, the finite temperature procedure is the same as in the previous section with the difference that the 
differential loops H and II have to be calculated at finite temperature. They now write 

E{U,U}(T,9 1 ,9 2 ,6 3 ) 



wli dej{vk ' 6) ^la-r 



xUr Je^Oa^Uu p (9 3 ,9 4 , 9) , (4.1) 



n{U u U 2 }(T, 0i,0 2 ,0 3 ) 

= A E /^A,0)[ /( - A) - /( ^ + r )]e(l ^ +qpJ ' A) 

xU 1 , lph (9 1 ,9,8 3 )U 2 ,i ph (9 4 ,9,e 2 ) . (4.2) 
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The function /(e) is the Fermi distribution at temperature T. The finite temperature version of the flow equations 



for the vertices z and for the correlation functions \ are again given by equations ( 3.20 ) and ( 3.21 ) but with modified 
D sc 



"=+, 



D? C {T,6)= ^A,0)tanh( — 

and D AF 



[/(M)-/(2H-^A)]9(|2H-^A|-A) 
2(,-M) 

AF 



(4.3) 



(4.4) 



We see that now both shells v = +,— contribute to Dp for A > \fi\, unlike in the zero temperature case where 
v = + contributions were forbidden by the Fermi distribution. In the discretized version of the flow equations, one 
calculates 5 and II from expressions fl3.24| ) and (^25| ), but with 



B pp (ii,i2,i;l,T) 



(2tt) 2 ^ 



(4.5) 



and 



2 



(2-) 2 ,i 



53 / d6J{vk,6) 



[/(^A) - /(a, +qp J] e (|^ +qp J - a) 

^ ASk„+q ph 



(4.6) 



The equations for finite temperature zf(T, and Xi(T, ii,t2) are ( 3.28 ) and ( 3.29 ) with -D; (T, i) calculated from 



( fU0| ), but using (Q and (gj). 

To find each point of the phase diagram, we have to find the fixed-point (A —> 0) value of the maximal eigenvalues 
xt(T,l) of the correlation functions xf(ii,«2); (<$ = SC,AF), It means that the complete renormalization from Z = 
to Z — > oo has to be done for each temperature. The flow of the quantities xf C (T,l) and xf F {T,l) is shown in 



figure 18 for several temperatures and three different values of the chemical potential. The susceptibilities for the 
non-interacting (U — 0) case are also shown. For all calculations the initial interaction was Uo — 4t/3 and we have 
cut the Brillouin zone into m, = 32 patches. The symmetry of the dominant superconducting correlations is for all 
cases B\ (which transforms as d x - y 2) and the dominant antiferromagnetic correlations have A\ (s) symmetry. They 
correspond to the symmetries of the strongest attractive components of the effective interactions y sc > AF found in 
the previous section. 

Let's concentrate first on figure |l8|(a) that shows the flow at half filling. The entire flow is in the parquet regime: 
the nesting is perfect. In the beginning of the flow where A; 3> T all correlation functions, bare or with correlations, 
scale as if the temperature was zero, i.e. like ln 2 (Ao/A) — I 2 . As the cutoff approaches the temperature, the flow 
starts to saturate. At the same time the effects of the interaction become more and more visible as we decrease the 
temperature. For all temperatures xf F an d Xc ° are enhanced from their bare values Xo F an d Xo C that are equal 
at half filling. As we approach the temperature T ss 0.0163 from above, the difference between the bare and the 
interacting cases increases rapidly; we interpret this temperature as the critical temperature. We have approximated 
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FIG. 18. The flow of the correlation functions with interaction (thick lines) and without interaction (thin lines) (a): at half 
filling for T/At = 0.03(A), 0.0204(B), 0.0163(C); (b): at \fi\/At = 0.002 for T/At = 0.0228(A), 0.0108(B), 0.0086(C); and (c): 
at \fj,\/4t = 0.006 for T/At = 0.03(A), 0.006(B), 0.0026(C). 
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FIG. 19. The temperature dependence of the fixed-point correlation functions for three different values of the chemical 
potential. 



the fixed point values Xc S (T) and Xo{T) of xi(T, I) and Xo(T, I) with their value for I = 10 (the corresponding energy 
A; is much smaller than any physical energy scale) . Figure |l9|(a) shows the temperature dependence of the fixed point 
values at half filling. The bare susceptibility scales as ln 2 (Ao/T). The interaction makes both susceptibilities diverge, 
but the antiferromagnetic one diverges first. 



Now we increase the chemical potential to |/u|/4£ = 0.002 (figure |1 8Kb)). The beginning of the flow, where A; 3> 
is still square-logarithmic. When I becomes close to — 6.9 the bare antiferromagnetic correlations start to be 
weaker because we approach the crossover from the parquet to the BCS regime. The non-analytic point is at I = 



at this point the flow equations ( p.28j ) and ( p.29| ) have a peak because of the van Hove singularity in Df(9) at 
9 = 0, 7r/2, 7r, 37r/2. Again, as we approach the critical temperature the effects of the interaction become stronger and 
stronger so that the difference between xt(T,l) and xt(T,l) increases more and more. The temperature dependence 
of the fixed point values for all correlation functions for the present case are shown in figure |l~9|(b). The critical 
temperature is T^ G /At = 0.0075, which is higher than |/x|/4f = 0.002. This means that the instability is still in the 
parquet regime, but not too deep: the proximity of the crossover already affects the antiferromagnetic correlations 
which start to loose their strength with respect to the superconducting correlations near the instability. However, both 
are still strongly enhanced and their flow is dominated by the parquet part (I < l^) for all temperatures T > T^ G . 

Let's increase further the chemical potential to \n\/4t = 0.006 (figure |l^(c)). The flow of the antiferromagnetic 
correlations saturates in the BCS regime (/ > Z p = 5.8), but both correlation functions SC and AF remain enhanced 
from their bare values. In the temperature dependence of their fixed point values (figure |l^(c)) one sees that only 
the superconducting instability is possible. The critical temperature is lower than the chemical potential, i.e. the 
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instability is in the BCS regime, in which the p-h part of the flow is negligible. The antiferromagnetic susceptibility 
even starts to decrease with the temperature when T <J (a. This happens because we do not adjust the wave vector of 
the SDW to the best nesting (incommensurate SDW) but we keep it for simplicity at (n, tt). Nothing drastic would 
happen even if we have taken small deviations of the best nesting wave vector from (71", 7r): the susceptibility would 
saturate as the temperature decreases because the differential p-h bubble decays in the BCS regime at any wave vector 
with positive power of A. 

The universal function f(\fi\/T®) (equation |3.3l| ), which determines the dependence of the critical temperature on 
the chemical potential, is practically the same as the one obtained in the zero temperature formalism. The final phase 
diagram is the one on figure At half filling the antiferromagnetic fluctuations are dominant over the supercon- 
ducting ones but both correlation functions diverge at the critical temperature. Upon doping the antiferromagnetic 
correlations loose their strength while the superconducting correlations remain strongly divergent. The divergence of 
AF correlation functions is completely suppressed if the critical temperature is in the BCS regime. 



V. CONCLUSION 

We have formulated the exact Kadanoff-Wilson-Polchinski (KWP) renormalization group for a general problem of 
interacting fermions on a two-dimensional lattice. In principle the generalization to higher dimensions is trivial. The 
procedure of the KWP renormalization scheme is to integrate out successively the degrees of freedom starting from 
high energies and to follow the renormalization of all terms in the effective action. We parametrize the renormalization 
by a high-energy cutoff A = Ao cxp(— /) determining the ring ±A around the Fermi surface. In order to take the 
whole Brillouin zone into account, the cutoff A is taken to be equal to the bandwidth (Ao = B.W.) at the beginning 
of the renormalization (1 = 0). As one proceeds with the mode elimination, vertices of all orders are created. To 
follow the exact renormalization of the effective action we need to know the flow (the dependence on I) of all vertices. 
The Polchinski equation (equations 2.16, 2.17 and figure ||) determines the differential flow of all vertices as functions 



of energy-momenta. In principle the fixed point solution (I — ► 00) of this equation gives us the exact connected 
Green functions of the modelES Clearly, the exact integration of Polchinski equation is impossible and for concrete 
calculations we have to truncate the effective action. 

The truncation at the sextic term (at the three-particle interaction term) generates the one-loop renormalization 



group for the two-body interaction. The truncated effective action is given by expression (2.21). Its renormalization 



is determined by the flow equation for the two-body interaction and for the self-energy. The flow equation for the 



interaction Ui is made of all one-loop diagrams bilinear in Ui as shown on figure ^ and by equation ( 2.25 ). Note 



that Ui is renormalized as a function of three energy-momenta (the fourth is conserved), i.e. this is a functional 



renormalization group. The (3 function (2.25) contains the contributions from the particle-particle (p-p) (/? PP ) and 
the particle-hole (p-h) diagrams (f3 p h)- The first term is called the BCS contribution in the literature, the next three 
terms are the zero sound (ZS) contribution, and the last term is the ZS' contribution to the differential flow. The flow 



equation for the interaction is not local in I as one can see from equations (2.28) and (2.29) for differential p-p and 



p-h bubbles: at some step I of the renormalization, Ui is renormalized by the values of U at former steps l pp (k, kj, IC2) 
and {^(k, ki, kg). This non-locality is the price we have to pay if we want to keep all contributions, logarithmic or 
not, that renormalize the interaction. In this way one takes correctly into account, for example, the p-h flow because 
of the imperfectly nested Fermi surface. The standard (local) Wilsonian RG0 that takes into account only dominant 
logarithmic diagrams (those with I = l pp = l p h) can give useful results only for the perfectly nested (but not square) 
Fermi surfaces or the Fermi surfaces far from being nested, so that the p-h part is negligible. 

We have applied the one-loop renormalization group to the Hubbard model on a square lattice near half filling. 
The interaction function U that we renormalize is dependent on the angular (Q) position of three momenta on the 
square Fermi surface (the fourth one is conserved). All radial momentum dependencies and energy (u>) dependencies 
are irrelevant to the Fermi liquid scaling. It is important that we allow variables 9 of the interacting particles to 
be anywhere on the square Fermi surface and not only in the configurations which give perfect nesting or zero total 
momentum (see figure ^). This means that we do not limit ourselves to the leading logarithmic parts of the flow but 
that we take all non-logarithmic contributions into account. 

From the explicit scale dependence of the differential flow for U we see two renormalization regimes (see figure [To|) . 
In the first regime, A; > We call it the parquet regime because both p-p and p-h contributions are important. 
The other regime exists in the non-half filled case when A; < |/z|. There, only p-p loops have a strong logarithmic 
flow while the p-h part decays to zero. We call this regime the BCS regime. The effective phase space (|£o(k)| < A) 
in the parquet regime is open so that the nesting is relevant (see figure ^|), while in the BCS regime the phase space is 
a closed regular ring of degrees of freedom around the Fermi surface so that perfect nesting is impossible (see figure 
The flow in the parquet regime is characterized by a strong coupling between the p-p and the p-h channels of 
renormalization. This coupling comes into play over the interactions that have a strong flow from both p-p and p-h 
diagrams. For the case of the (nearly) square Fermi surface these are all interactions between electrons from opposite 
sheets of the Fermi surface. 

The leading correlations in the Hubbard model are expected to be antiferromagnetic and/or superconducting. To 
give a precise answer to this question, we use Kadanoff-Wilson-Polchinski procedure to construct the renormalization- 
group equations for the angle resolved correlation functions Xi (^1^2) and Xi (^1^2) for antiferromagnetism and 



superconductivity, defined by equation (3.10). At a given step A(Z) of the renormalization these correlation functions 
measure the linear response of the electrons outside the shell ±A around the Fermi surface. We take the static 



long-wave limit. The renormalization equation for xf F ' SC is ( 3.21 ). The renormalization of the correlation functions 



depends on the renormalization of the vertices z AF ' SC (eq. 3.20). Furthermore, from equations ( 3.2C ) and (3.21) 
one sees that the flows of the susceptibilities and of the vertices depend on the flows of the corresponding effective 



interactions Vf c and V^ F , given by (3.5) and ( 3.16 ) respectively. 

The flow equations for the interaction Ui (2.25|), for the vertices z x ' (3. 20]) and for the correlation functions 



Xi ' F (3-21) can be integrated numerically if we discretize their ^-dependence. The coupling function is then 
approximated by a set of coupling constants. The vertices and correlation functions become discrete matrices. Using 
physical and geometrical symmetries we reduce number of the coupling constants to a set of the independent ones (see 
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figure [lj). The functional renormalization-group equations become a set of equations, one for each coupling constant 
and for each matrix element of zf C,AF and xf C ' AF . 

We have solved the renormalization equations for up to m.; = 32 angular patches. The typical flow of the coupling 
constants is shown on figure [l5|. There is a critical scale for which a large number of coupling constants diverge. We 
associate the corresponding cutoff to the critical temperature T^ G ■ Its dependence on the chemical potential is shown 
on figure |l6| (solid line) together with the RPA result (dashed line). At the line T c flG (/x) the electronic correlations 
are strongly enhanced. The type and the form of the corresponding microscopical fluctuating fields are given by the 
dominating eigenvalues (and their eigenvectors) of the correlation matrices xf C ' AF . These are determined by the 
dominant attractive eigenvalues V^ c and V AF of the effective interactions. For all values of the chemical potential 
studied in this work, the eigenvalue Vf c corresponds to d x 2_ y 2 (or B\) singlet superconductivity while V AF is an 
s-wave ( Ai representation) . The flow of the interactions Vf c and V AF in the vicinity of the critical point A = T c fiG 
is shown on figure [l7]. At half filling V AF is dominant. Upon doping, the divergence of V AF loses its strength and 
the divergence of V c becomes dominant. 

To determine more precisely the dominant fluctuations near T c flG , we have done one further step in the 
renormalization-group formalism: we have introduced the temperature explicitly into the flow equations. In this 
formalism the cutoff A has no more a physical meaning of the effective temperature. At given temperature T the 
physical information is contained in the fixed point (A — > 0) of the correlation functions. This extension of the 
formalism was necessary because in the zero-temperature procedure it was not possible to have the divergence of 
xf F at the same scale l c as the divergence of the coupling V AF : for any nonzero chemical potential the flow of 
xf F has a finite retardation in I (see equations ( 3.20| ) and ( 3.21 ). Thus xf F diverges later, aX I > l c . In the finite 
temperature formalism x AF \ xf G > ^/ and Vf c all diverge at the same temperature. The price to pay is that for 
each temperature we have to integrate the complete flow all the way from I — to I = oo and to follow how the result 
changes with the temperature. 

The flow of the dominant eigenvalues of x AF and xf C at a few different values of the temperature is shown on figure 
|l8| for three different values of the chemical potential. Both correlation functions are always enhanced with respect 
to their bare (U = 0) values. The temperature dependence of the fixed point correlation functions is shown on figure 
|l9| . The critical temperature T C RG found by the finite temperature method is practically the same as the one found 
by the zero-temperature calculations, but now we are able to follow explicitly the enhancement of the correlations 
of both types in the vicinity of the instability. It is clearly visible how the doping favors superconductivity and how 
the divergence of x* AF is completely suppressed if the instability is in the BCS regime, i.e. if T^ G < This result 
justifies the phase diagram on figure [D| 

In the low-doping regime, both correlations are strongly enhanced and the low-temperature phase can be in 
principle a mixture of both (quasi-)long-range-orders, with the superconducting component falling to zero at half 
filling. The instability is in the parquet regime: the critical fluctuations are a mixture of two fluctuating channels and 
can not be treated by an effective mean-field theory like BCS or RPA. In other words, the parquet regime is deeply 
non Migdalian so that the vertex corrections are as important as the p-p loops. The vertex corrections can not be 
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seen as small any more but have to be taken at all orders, just as the p-p diagrams, and together with other p-h loops. 

The situation is less complicated in the BCS regime. There, a low energy effective action can be constructed so that 
only the p-p diagrams contribute to the instability while the p-h parts (antifcrromagnetic tendencies) are irrelevant. 
The attractive d-wave component of the Cooper amplitude in the LEE A is due to the higher energies (e > \fi\) where 
the p-h diagrams are important. In the BCS regime only superconductivity is possible as a low-temperature order. 

As we are considering a two-dimensional system, one should be careful about the interpretation of T c : in the case 
of magnetism, this indicates the onset of well-defined finite-range correlations. For weak interactions, this is typically 
a very well-defined crossover .0 In the case of pairing T^ G can be identified with the onset of quasi-long-range order. 
Furthermore, the line between AF(SCd) and SC phases on figure [l^ is only partially determined in our calculations: 
we only know that at temperature near , this line is close to the crossover line T = but at lower temperatures 
we can not say anything about its position. 

It is difficult to discuss the experimental results from the point of view of our phase diagram. First of all, the 
one-loop renormalization-group is a weak coupling perturbative method while the interactions in the copper-oxide 
superconductors are moderate-to-strong. For that reason our phase diagram can be compared to the experiments 
only qualitatively. Furthermore, in our calculations we have neglected self-energy corrections, which are in Polchin- 



ski's formalism given by Hartree-Fock-like terms with renormalized to- and q-dependent vertices (eq. 2.30). The 
broadening and redistribution of the spectral weight of the quasiparticles is then determined by the dynamics of the 
vertex, which is irrelevant and is therefore neglected. One should however notice that at the two-loop level self-energy 
effects become important, as known from the one-dimensional case.0 In that sense, our T c flG should be understood 
as a temperature where the effects of interactions start to change strongly not only the two-particles correlations, 
but the single particle properties as well. For that reason it seems natural to associate the temperature T^ G to 
the crossover temperature T co found in the cuprates. The parquet regime would then correspond roughly to the 
under-doped situation and the BCS regime to the over-doped regime. 

The "phase" AF(SCd) corresponds then to the antifcrromagnetism and to the pseudo-gap regime: the antiferro- 
magnetic correlations and the localization tendencies are there accompanied more and more with the superconducting 
correlations as we approach the crossover line T = |/z|. We expect that the critical temperature for antiferromagnetism 
and for superconductivity in this regime is lower than T^ G because of the self-energy corrections. In other words, at 
temperature T^ G in the parquet regime, the local antiferromagnetic moments and d-wave singlets are created with 
finite correlation lengths. This gives rise to the pseudogap in both spin and charge responses, together with the pre- 
cursors of both antiferromagnetism and d-wave superconductivity. The long-range-order between superconductivity 
and antiferromagnetism is perhaps absent due to the fact that both types of fluctuations are strong. That is the 
central idea of the SO (5) modelsEI for the high-T c superconductivity. In that language our T^ G would play the role 
of the mean-field critical temperature. 

In the BCS regime only the superconducting fluctuations are critical. We associate thus the phase SCd to the 
overdoped regime. From large N argumentsE3@ we know that the self-energy corrections disappear as T c /t if the Fermi 
surface is not nested. This is the case in the BCS regime where the nesting processes are irrelevant. Consequently, 
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the critical temperature in this regime is well approximated by T^. This is in agreement with the experiments: in 
the overdoped regime the crossover temperature T co is equal to the critical temperature for the superconductivity. 
Finally, mean-field argumentsil suggest that one expects an incommensurate SDW (ICSDW) only in the BCS regime 
and only where the imperfect nesting is still strong, i.e. not far from the crossover T — However, the precision 
of our calculation (we cut the Brillouin zone into up to 32 ^-patches) is not sufficient to check whether a magnetic 
correlation function diverges at some incommensurate wave vector. In any case, the incommensurate SDW and d- 
superconductivity are not in competition because they appear at different places on the Fermi surface: SCd in the 
corners and ICSDW on the flat parts; one thus expects their coexistence. 

Altogether, the phase diagram on figure |l6| has important similarities to the experimental phase diagrams. The one- 
loop renormalization-group, taking into account electron-electron and electron-hole processes on the same footing 
reveals the essence of the physics of a doped half filled band of correlated electrons. 
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APPENDIX A: INTERACTION (7(1,2,3) AND ITS SYMMETRIES 

The most general spin-rotation invariant interaction term can be written in several ways. One way is in terms of 
charge-charge and spin-spin interactions 

V C {K X ,K 2 ,K 3 )C{K 2 , Ki)C{K 3 , K ± ) + U a {K u K 2 , K 3 )S{K 2 , K 4 ) ■ S(K 3 , K x ) , (Al) 

where C et Si are 

C{K 3 ,K x ) = Y,*°k^°Ki ! S i (K 3 ,K 1 ) = J2*<rKA*'*«'Ki ■ (A2) 

The summation over all three energy-momenta [K\, K2, K 3 ) is assumed and K± = K\ + K2 — K 3 . On the other hand, 
the interaction can also be written as a sum of one term with equal (ct = ct') and one with opposite (ct = — ct') spin 
quantum numbers, with corresponding coupling functions named U» (K\, K 2 , K 3 ) and U±(K\, K2, K 3 ): 

U ll (K ll K2,K 3 )^ aKa ^ aKi ^ aK . 2 ^ aKl +U A _{K 1 ,K 2 ,K 3 )^ (TK3 ^_ aKi ^_ a . K2 ^ lTKl (A3) 

with the summation over spin indices assumed. Spin-rotation invariance allows us to write the interaction part of the 
action as a sum of the singlet (|c? + a'\ = 0) and triplet (|ct + ct'| = ^/2) parts: 

s(K 4 , K 3 )U s (K ly K 2 , K^siK^Kx) + i^(K 4 , K 3 )U A (K 1 ,K 2 , K 3 )t^K 2 , Ki), (A4) 
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where s and are the variables of annihilation of the singlet and triplet states 
s(K 2 ,Kx) = <r* ff jc a 

* a 



(A5) 



(A6) 



All coupling functions in the equations (A2,A3,A4) possess the symmetry related to momentum-exchange and time- 
inversion. If T(K±, K 2 , K 3 , K 4 ) is a coupling function two exchange operators can be defined as 



XT{K X ,K 2 , K 3 ,K 4 ) = T{K 2 ,K X ,K 3 , K 4 ) 



and 



X^(K 1 ,K 2 ,K 3 ,K 4 )=^(K 1 ,K 2 ,K 4 ,K 3 ). 
The time inversion operator T is 

TT(K U K 2 , K 3 , K 4 ) = T(K 3 , K 4 , K U K 2 ). 
The symmetries of the coupling function are the time inversion symmetry 

TT = T 
and the exchange symmetry 

XXT = T . 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



Both symmetries can be easily checked for the coupling functions in expression (A3). We will see that all other 
couplings can be derived from U± only and have the same symmetry properties upon X and T operations. It is easy 
to see that XT = XT if TT — T: exchanging particles 1 and 2 or particles 3 and 4 are exchanged is equivalent. 

We want now to find the relations between the six coupling functions in equations (A.2, A3|,A4). Using the Pauli 
principle one gets 



U\\ =u c + u ai 

U^ = U c -U a - 2XU a . 
Let's suppose U c and U a to be two independent functions. We can write them in the form 
U C =^(2-X)U 1 + U 2 , 



X 



(A12) 
(A13) 

(A14) 
(A15) 



If we now choose U\ = U± it follows from (A13) that U 2 = 0. It means that the most general interaction can be 
written in terms of a single function U±, without losing generality. The function J7|| is also contained in U±. Namely, 
from two equal-spin electrons one can build only a triplet state (antisymmetric under X) so that 

jA 



[7|l = U 1 



(A16) 
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while 

U±_ = U A + II s , (A17) 

containing the singlet and the triplet interactions. U A and U s can be seen as the antisymmetric and symmetric parts 
of the same function. This function is simply U± . 

We see that all coupling functions are contained in U± which we call simply U or Ui to make its scale dependence 
explicit. One thus have 

U C = ^(2-X)U, Ua = -^U, (A18) 

U A = U {] =±(1-X)U, U S = ±(1+X)U. (A19) 

The effective coupling function for the renormalization of the AF correlation function (eq. 3.16| ) is obtained from 
the spin coupling 

V AF {9 1 ,9 2 )=W al {9 u 9 2 ,9 1 ) , (A20) 

where we take only the ^-dependence of the coupling functions into account. The angle 9 is related to the angle 9 in 
such a way that the momentum difference between the particles k((9) and k(#) is the perfect nesting vector (±7r, ±7r). 
The coupling function for the CDW at q = (tt, it) would be 

V l CDW (9 1 ,9 2 )=4U cl (9 1 ,9 2 ,9~ 1 ) . (A21) 
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